Construction of multiple spherical branes cosmological scenario 
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The general solution of Einstein's gravity equation in D dimensions for an anisotropic and spher- 
ically symmetric matter distribution is calculated in a bulk with position dependent cosmological 
constant. Results for n concentric (D — 2)— branes with arbitrary mass, radius, and pressure with 
different cosmological constant between branes are found. It is shown how the different cosmolog- 
ical constants contribute to the effective mass of each brane. It is also shown that the equation of 
state for each brane influences the dynamics of branes, which can be divided into eras according to 
the dominant matter. This scenario can be used to model the universe in the D — 5 case, which 
may presents a phenomenology richer than the current models. The evolution law of the branes 
is studied, and the anisotropic pressure that removes divergences is found. The Randall-Sundrum 
metric in an outside the region in the flat branes limit is also derived. 
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I. INTRODUCTION 



The general model for the cosmos is based on the description of the universe as a perfect fluid that admits a global 
cosmic time. This scenario has a space-time with constant curvature given by the Fricdmann- Robertson- Walker metric, 
where its dynamics is determined by a cosmological scale factor that depends on the fluid state equation. Despite the 
success of that model in the description of the primordial nucleosintesys and the cosmic microwave background, it 
has failures that led to the emergence of new models. Among the major flaws of the current model are the problem of 
dark energy showing the accelerated expansion in the currently observed universe, and the dark matter which is the 
divergence between the rotation of the halo of some galaxies and the amount of matter contained in them according 
to gravitational dynamics [l| . 

Cosmological models with extra dimensions appeared first in Kaluza-Klein models with extra dimensions and 
later in Randall- Sundrum scenarios [H, Q . These models describe the observed universe as a brane universe in a 
hyper-dimensional space-time. Despite the fact that the Friedmann-Robcrtson- Walker metric did not determine the 
geometry of the observed universe. The majority of studies focused on plane geometry. Because of its simplicity, this 
geometry is not able to change the dynamics of the universe and thus cannot solve the problem of dark matter or the 
initial singularity. 

Although the first studies to describe the universe as a spherical shell back to the 80 's 041], the spherical brane- 
universe has shown very rich phenomenology in the past decade 0, @] ■ Besides being compatible with the observational 
data [9l4Tl|. the models provide an explanation for; the galaxy isotropic runaway (isotropic expansion), the existence 
of a preferred frame, and a cosmic time. They show how the introduction of different cosmological constants in 
each region of the bulk can change the dynamics of the cosmological scale factor so as to make it compatible with 
the observed dynamics fl2l [l3j without the introduction of dark energy [lij ■ Similar to other models with extra 
dimensions, the spherical shell models open the possible to obtain an energy scale in order to solve the pro blem of 
hierarchy |15| and can be used as a basis for systems with varying speed of light in the observed universe |16j . 

The introduction of other branes and different cosmological constants can modify the overall dynamics of the 
observed universe. Local density fluctuations of density can change the local dynamics such as galactic dynamics 
(since the field of other branes interacts gravitationally with the matter of the brane-universe) without dark matter. 
Herein this piece of work we extend and generalize the scenario of the world as one expanding shell to multiple 
concentric spherical (D — 2)-branes in a D dimensional space-time. For this, we solve the Einstein's equation in D 
dimensions to n (D — 2) —branes with different masses in a space with different cosmological constants between the 
branes. A previous study considered a continuous distributions of matter. However, only one cosmological constant 
was used We solve the I?— dimensional case, but for a cosmological model we limited ourselves to the case D = 5, 
since the observed universe has only three spatial dimensions. 

This work is organized as follows: In the second section we review the Einstein's equations in D dimensions 
with a cosmological constant for spherically symmetric matter distribution. In the third section we solve this set 
of equations for n shells with different cosmological constants A between them. In Sec. 4, the energy-momentum 
tensor conservation law is used to determine the possible anisotropic pressure which removes the divergences in brane 
evolution equation. In the fifth section we particularize the solution found to take the flat brane limit in order to 
obtain the Randall-Sundrum metric in the exterior region. In the last section we discuss the conclusions and possible 
consequences. 



II. STATIC AND SPHERICALLY SYMMETRIC SPACE-TIME IN D DIMENSIONS 



To learn about the gravitational effect of a distribution of matter we must determine the geometry of space-time. 
For this we need to know the D(D + l)/2 independent components of the metric solving the Einstein's equation. 
However, it is possible to use the symmetry of the problem to reduce these components to just two, given by the 
invariant line element Q, 

ds 2 = -A(r, t)dt 2 + B(r, t)dr 2 + r 2 di} 2 D _ 2 (1) 

where is the element of solid angle in D dimensions, formed by D — 2 angular variables. 

Therefore we are left only with two functions, A(r, t) and B(r, t), to be determined by the Einstein's equation in D 
dimensions 

r;-Irsz + a5; = k d t;, (2) 

where A is the cosmological constant, which depends on r and possibly on t. Also kjj is the gravitational coupling 
constant in D dimensions. Due to the symmetries of the problem we only have four non null independent components 
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of the Einstein's equation (|2|) , which are 
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where the prime means derivation with respect to r and the dot is the derivative with respect to t. 

We can see that if we know Tq, and A we can, from ((3]) and completely determine the solutions with two 
boundary conditions. This comes from the fact that we have two first order differential equations. In this case the 
remaining equations determine the flow of energy Tq 1 , and the tangential stresses To find the exact solution we 
need to specify the form of matter Tjf which we use. 



III. GENERAL SOLUTION FOR THIN SPHERICAL BRANES 



The cosmological scenario we shall consider consists of n concentric spherical delta type (D — 2)-brancs in a D di- 
mensional space with different cosmological constant between them. As said in the introduction this is a generalization 
of • For this we fix the energy-momentum tensor and the cosmological constant to the form 

n n 

7 1 o (^^) = -E^( r -^)' Tt = Y,Pi5(r-Ri), (7) 

i=l i=l 

and 

A(r, t) = {D " 1} 2 (g ~ 2) ]T Xi [0(r - R.) - 9{r - R l+1 )] , (8) 

i=0 

where the dependence on t should be solely due to the branes radii (i?.; = Ri(t)). The 9 function is defined in such 
way that is 1 when the argument vanish, it's made to ensure that above expresion cover all space, including r = 
point. 

The cosmological constant can be understood as a special fluid, so we can think that the difference between the 
cosmological constant is because each brane contains a fluid with different density. Fixing Tq and A we can find 
B(r,t) using equation © in the form 

«d'$ = -§^[t»^(1-B- 1 )] , + A (9) 

and according to the above equation, B has a first order discontinuity in Ri because Tq has a second order one and 
A has first order discontinuity only. Where we consider that in the region Ri < r < B(r,t) = Sj(r), since in 

this region B does not depend on t. The above-mentioned time dependence occurs in the region where this solution 
is valid. 

The region between the branes has no matter. Therefore, the equation (|5|) assures us that the solution is static in 
this region. This information is contained in the Birkhoff 's theorem. Integrating ((9]) from Rj — e to Rj + e and taking 
the limit e — > 0, we obtain the discontinuity in the point r = Rj 

Bj\ Rj ) Bj\( Rj ) = -j^PjRi (10) 



The limit e — > eliminate the A term because its divergence is first order only. Integrating @ from Rj +e to r < Rj+i, 
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where B is continuous, and taking the limit e — > we obtain 
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By recurrence we find that 
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where AAi = A, — Ai_i. Considering that inside all branes the solution is a de Sitter vacuum, i.e., Bq(t) 
(l — A r 2 ) , we get 
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The above solution is valid only in the region Rj < r < Rj+i, but we can write the solution valid in any region in the 
form 



B- l (r,t) = l- 
where M(r, t) and A(r, t) are denned by 
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and the time dependence is implicit in Ri. It is important to note that M(r,t) is not positive defined, in order to 
enable a repulsive gravitational situation. Using the above definition in ([¥]) we find the equation which governs A 
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Taking the way T^ 1 was fixed at ([7]), A has a second order discontinuity. Now, using the same procedure to find B we 
can show that 



where 



A j (r) = Bj\r)A {R{)B a (R 1 ) f[ ^^V ' ■ 

i=X ^ ' 



7r» = ^^RiBiiRAPi. 



The asymptotic behavior of B(r,t) is lim^oo B(r) = [l — A(r)r 2 ] , which is the generalization of the de Sitter 
vacuum to a cosmological constant that is position dependent. Likewise we expect that A(r) behaves asymptotically 
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as the vacuum, i.e., lira r _^oo A(r) = 1 — A(r)r 2 , so that we can use this to fix the multiplicative constants appearing 
in the temporal solution and write 

In the same way we did for we can rewrite the above solution in order to be valid in all space as 

n 

A{r,t) = B- x (r,t) JJ e -»«o(n«-r) x 

JRA \ 1 

(13) 



e 



where Bj is defined by (fTTj) . 

The solutions (fT3")) and (fT2")) are generalizations of the Kottler solution [l8| in £) dimensions with position dependent 



A. In the case where A is constant, these solutions agree with those found by Das [17|. However, these solutions only 
make sense if the branes are not in a time-like region. In order to avoid a singularity in the solutions, we impose that 
the radius of the brane relates to the masses so that they are beyond their respective generalized Kottler's radii, i.e. 

- XiR^Rf- 1 + Rf-z - 2GM(Ri) > 0. (14) 

The above solutions perfectly agree with the Birkhoff's theorem, and despite a constant, it is the Schwarzschild 
solutions with a cosmological constant (Kottler Solution). The temporal dependence of the solutions is in i?,;, so that 
in each region the solution is static. The multiplicative constants in the temporal part of the solutions indicates the 
gravitational redshift, even within the shells. Mathematically this makes the solution continuous in all regions. 

The time dependence on B is given exclusively by Therefore in a dynamical case (Ri = Ri(t)) we can obtain 
the energy flow from ([5]). It's easy to show that 



T^-^piVidir-Ri), 



i=i 



where Vi = R4. The tangential stresses can be obtained from (|6]), but it is easier to compute from the energy- 
momentum tensor conservation law. 



IV. ENERGY-MOMENTUM TENSOR CONSERVATION LAW 



The Einstein's field equation relates the energy-momentum tensor and the metric tensor. But due the symmetry 
only two components of energy-momentum tensor are necessary to determine the metric. So the other terms of energy- 
momentum tensor are determined by Einstein's equation or by a conservation law. Taking the covariant derivative in 
the equation ([2]) we obtain the D dimensional conservation law. 
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(15) 



the above equation states that the energy and momentum are not conserved inside the brane because the extra 
dimensional pressure given by the difference between cosmological constants. This difference can be used to model 
the dark energy, which makes the universe expand. But in our case no strange matter in needed inside the brane like 
the usual dark matter models. In terms of independent components the above conservation law can be written as 
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The first equation is trivially satisfied if we use the known solution (fT3|) and (fT2)) and energy-momentum tensor 
components. The second equation gives us the propagation speed of the brane as a function of the tangential stress, 
the masses, and the cosmological constant in the same way as in (j6|). Taking 

n 
4=1 

we can integrate the unsolved component of the conservation law from Rj — e to i?j + e to obtain 

1 [piVi + pM] + — — [Pi - T t ] + 
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as the functions A and B have second order divergences in r = Ri the last term in above expression have the same 
divergence. Analyzing separately the divergent terms 
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To avoid a real divergence we need to fix 
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The first case indicates a cosmological constant state equation type, this equation is the only state that is independent 
of motion, i.e. the properties of a fluid with this state equation is independent of its movement. Therefore, it was 
already expected that the divergences found in the dynamic case can be removed. The second case relates the normal 
pressure with the brane velocity. That indicates an increase in the pressure if the velocity increases to keep the 
spherical shape of the brane. This relationship ensures that P could vanishes in the static case, as in the Randall- 
Sundrum scenario. Assuming a linear state equation relating the tangential stresses and the energy density, Tj = jiPi, 
and defining, for the i-th brane, the time 



dU = v / A l {R l )/B l {R l )dt 



the brane evolution is given by 
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where Ui = This indicates a different dynamic for each cosmological eras driven by a different state equation , 
i.e., by the related tangential pressures. 
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V. THE RANDALL-SUNDRUM FLAT BRANE LIMIT 



In the previous sections we found the general solution to n spherical branes in a £>-dimensional space-time with 
different cosmological constant between them. To find a scenario similar to Randall-Sundrum we need to fix D = 5, 
n = 1 and Aq = A x . In this case the exterior solution is 



ds 2 



-f(r)dt 2 + f-\r)dr 2 + r 2 dSl 2 , 



where 
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In order to obtain the Randall-Sundrum metric we define 
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where A = — k 2 to avoid the de Sitter horizon. To obtain the fiat brane limit we will consider that R and consequently 
r tends to infinity. In this limit the dominant term, regarding that M grows as R 3 , is 



Z= 2l ln 



R 2 



writing r as function of z, we obtain the line element 



ds 2 = -k 2 e 2kz R 2 dt 2 + dz 2 + e 2kz R 2 dn 2 . 



Putting the constants into coordinates we obtain the Randall-Sundrum metric 



ds 2 = 



2k; 



r\ >lv dx p 'dx v + dz 2 



The exponential in the warp factor is positive because the bulk is anti de Sitter, instead the original RS scenario. 



VI. CONCLUSIONS AND PERSPECTIVES 



In this work we built a scenario of multiple concentric membranes through the solution of Einstein's equation in D 
dimensions with different cosmological constant in each region. The results we found may serve as a basis for more 
specific scenarios, through the fixation of radii, masses on each brane, and bulk cosmological constant. In a dynamical 
case, the solutions we found can be used to model the universe for D = 5. The model used here is more accurate 
than the previous ones, and as a consequence we have a multiplicative constant that appears in the temporal solution, 
which is the redshift measured by observers in the region inside the brane. 

Through the momentum-energy tensor conservation law we obtained two possible anisotropic pressures that remove 
the singularities in the branes dynamics. These two possible pressures give us two possible fixations leading to more 
freedom in the construction of a cosmological model that can better fits the observed data. The tangential pressure 
and the difference between the cosmological constants are responsible for the evolution of each brane according to 
Eq. (|18|) . This tangential pressure is found from a state equation that is determined by the dominant matter in 
each cosmological era. We show that the difference between the cosmological constants modifies the effective mass of 
the matter distribution, and can be fixed in a way that the observed universe expansion rate is independent of the 
dark energy. Finally, we were able to arrive at the Randall-Sundrum metric to a anti-de Sitter bulk from calculating 
the external solution in the limit of plane branes. This metric was first introduced in the literature as an ansatzQ. 
Now, it is derived from the Kottler anti-de Sitter solution. The main extension of the model developed here is the 
phenomcnological study of cosmology generated by solving a dynamic Universe Brane equation. 
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